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I .  Introduction 

This  report  contains  a  target  position  estimation  program  for 
each  of  the  following  calculators:  The  Hewlett-Packard  HP-4 lev, 
the  Sharp  PC-1500  (or  Radio  Shack  TRS-80  PC-2) ,  the  Sharp  PC-1211 
(or  Radio  Shack  TRS-80  PC-1) ,  the  Casio  FX-702P  and  the  Texas 
Instruments  TI-59. 

The  programs  provide  a  means  of  implementing  a  position  esti¬ 
mation  procedure  that  is  described  in  Appendix  1.  The  procedure 
is  based  on  the  following  assumptions:  Bearings  taken  on  or  from 
a  target  are  available  from  two  or  more  stations  of  known  position. 
The  positions  of  the  stations  and  the  target  are  such  that  they 
can  be  considered  to  be  on  the  surface  of  a  plane  (a  flat  earth 
model) .  The  error  in  a  bearing  taken  on  or  from  a  station  is  a 
normal  random  variable.  Its  standard  deviation  (bearing  error) 
is  known  and  its  mean  (bias)  is  zero  (if  there  is  bias  in  a  bear¬ 
ing,  it  is  known  and  it  is  removed) .  The  errors  in  bearing  mea¬ 
surements  are  independent. 

The  position  estimation  procedure  requires  an  initial  estimate 
of  the  target's  position.  In  the  programs  that  are  given  here, 
the  initial  estimate  is  at  the  intersection  of  the  bearing  lines 
determined  by  the  first  two  target  bearings  that  are  input  to 
the  program.  Because  of  the  use  of  this  method  to  choose  an  ini¬ 
tial  estimate,  the  first  two  stations,  with  respect  to  the  order 
of  data  input,  should  be  chosen  so  that  the  intersection  of  their 
bearing  lines  is  likely  to  be  closer  to  the  target’s  position 
than  the  intersection  of  the  bearing  lines  from  any  other  pair 
of  stations. 
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The  target  ranges  and  errors  of  the  bearings  of  the  first  two 
stations  determine  the  distance  between  the  intersection  of  their 
bearing  lines  and  the  target's  position.  In  particular,  if  the 
angular  separation  between  two  stations  as  seen  from  the  target 
is  small  relative  to  the  bearing  error  of  one  or  both  of  the  sta¬ 
tions,  the  bearing  lines  from  the  two  stations  may  not  intersect. 
If  they  do  not  intersect  and  they  are  not  parallel,  their  reci¬ 
procal  bearing  lines  will.  In  this  case,  if  the  observed  bearings 
from  the  two  stations  were  the  first  two  bearing  inputs  to  the 
program,  the  initial  estimate  of  the  target's  position  would  be 
at  the  intersection  of  the  two  reciprocal  bearing  lines,  and  a 
gross  error  in  the  final  position  estimate  could  result. 


II.  Program  User  Instructions 

Before  using  a  program  for  the  first  time,  refer  to  the  notes 
that  follow  the  user's  instructions.  Also,  note  the  comments  in 
Section  I  regarding  the  relationship  between  the  accuracy  of  the 
computed  estimates  and  the  order  of  data  entry. 

For  a  PC-1211,  use  the  PC-1500  User  Instructions.  Use  the 
DEF  mode  and  substitute  SHFT  for  DEF  wherever  DBF  appears  in  the 
instructions. 

Station  positions  are  determined  with  respect  to  a  known  ref¬ 
erence  point.  The  reference  point  can  be  a  station  position,  in 
this  case  the  station  bearing  and  range  from  the  reference  point 
are  both  zero. 

Program  listings  are  given  in  Section  IV  and  two  example 
applications  are  discussed  in  Section  III. 
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Position  Estimation  Program 
HP-41CV  User  Instructions 


Step 

Instruction 

Prompt 

Press 

1 

Select  USER  mode  and  enter  the  program 
(see  Note  1) .  Press  £+  to  run  the 
program. 

E+ 

2 

Key  in,  in  decimal  degrees,  the  obser¬ 
ved  bearing  of  the  target  from  a 
station  or  the  reciprocal  of  the  ob¬ 
served  bearing  of  the  station  from 
the  target  (see  Note  2) . 

OBS  BR6? 

R/S 

3 

Key  in  the  station  bearing  in  decimal 
degrees  from  the  reference  point.  Use 
zero  if  the  station  is  the  reference 
point. 

STA  BRG? 

R/S 

4 

Key  in  the  station  range  (in  any  units) 
from  the  reference  point.  Use  zero 
if  the  station  is  the  reference  point. 

STA  RNG? 

R/S 

5 

Key  in  the  bearing  error  (standard 
deviation)  in  decimal  degrees. 

R/S 

6 

Repeat  Steps  2,  3,  4  &  5  for  one  or 
more  additional  stations. 

7 

Compute  bearing  &  range  estimates. 

1/x 

8 

A  computed  target  bearing  estimate  in 
decimal  degrees  is  displayed. 

BRG=est. 

R/S 

9 

A  computed  target  range  estimate  in 
meters.  For  an  elliptical  contain¬ 
ment  region,  go  to  Step  10  or  Step  17. 

To  enter  additional  data  from  new  or 
old  stations,  go  to  Step  24. 

RNG*est. 

10 

For  a  containment  ellipse  of  a  given 
containment  probability,  press  LN. 

LN 

11 

Key  in  the  desired  containment  proba¬ 
bility. 

PRB? 

R/S 

12 

Computed  value  of  the  ellipse  size: 

(See  Note  3.) 

SIZE*val. 

R/S 

Step 

Instruction 

Prompt 

Press 

13 

Computed  semi-major  axis  length: 

SMJ=val . 

R/S 

14 

Computed  major  axis  direction: 

DRC=val. 

R/S 

15 

Computed  semi -minor  axis  length: 

SMI=val . 

R/S 

16 

Computed  containment  ellipse  area: 

(See  Note  4 . ) 

A  =  val. 

17 

For  a  containment  ellipse  of  a  given 
size,  press  LOG. 

LOG 

18 

Key  in  the  desired  containment  ellipse 
size  (see  Note  3) . 

SIZE? 

R/S 

19 

Computed  value  of  the  containment 
probability: 

PRB=val . 

R/S 

20 

Computed  semi-major  axis  length: 

SMJ=val. 

R/S 

21 

Computed  major  axis  direction: 

DPC=val . 

R/S 

22 

Computed  semi-minor  axis  length: 

SMI=val . 

R/S 

23 

Computed  containment  ellipse  area: 

(See  Note  4.) 

A  =  val. 

24 

To  enter  additional  data  from  new  or 
old  stations,  press  /x  and  then  repeat 
Steps  2,  3,  4  &  5 

* 
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Notes : 


1.  The  program  size  is  35.  The  key  assignments  are:  TPE  -*■  E+, 

CON  ■*  /x,  EST  -*•  1/x,  SIZ  -*■  LOG  and  PRB  •*  LN.  If  they  are  not 
present,  they  must  be  made  in  order  to  follow  the  user 
instructions.  For  an  alternative,  see  Appendix  2. 

2.  Reciprocal  bearings  are  used  when  bearings  are  taken  on  known 
positions  from  the  unknown  position  (target) . 

3.  In  the  model  that  is  the  basis  for  the  program,  for  a  given 
probability  of  containment,  the  minimum  area  containment  region 
is  an  ellipse  centered  on  the  position  estimate.  The  semi¬ 
major  axis  =  kaMj  an<*  semi-minor  axis  =  koMI  where  k  is  the 
size  of  the  ellipse  and  and  aMJ.  are  the  standard  deviations 
(uncertainty  measure)  of  the  position  estimate  in  the  major 
axis  and  minor  axis  directions. 

4.  The  area  units  are  the  range  units  squared. 


Position  Estimation  Program 
PC-1500  User  Instructions 


Step 


Instruction 


Prompt 


Press 


STA  RNG? 


BRG  ERR? 


Enter  the  program.  To  run  the  program, 
press  DEF,  A.  For  a  PC-1211,  see  Note  5. 

Key  in,  in  decimal  degrees,  the  observed  OBS  BRG? 

bearing  of  the  target  from  a  station  or 
the  reciprocal  of  the  observed  bearing 
of  the  station  from  the  target  (see 
Note  1) . 

Key  in  the  station  bearing  in  decimal  STA  BRG? 

degrees  from  the  reference  point.  Use 
zero  if  the  station  is  the  reference 
point. 

Key  in  the  station  range  (in  any  units)  STA  RNG? 

from  the  reference  point.  Use  zero  if 
the  station  is  the  reference  point. 

Key  in  the  bearing  error  (standard 
deviation)  in  decimal  degrees. 

Repeat  Steps  2,  3,  4  &  5  for  one  or 
more  additional  stations. 

Compute  bearing  &  range  estimates. 

A  computed  target  bearing  estimate  in 
decimal  degrees  is  displayed. 

A  computed  target  range  estimate  is  dis-  RNG=est. 

played.  For  an  elliptical  contain¬ 
ment  region,  go  to  Step  10  or  Step  17. 

To  enter  additional  data  from  new  or 
old  stations,  go  to  Step  24. 

For  a  containment  ellipse  of  a  given 
containment  probability,  press  DEF,  X. 

Key  in  the  desired  containment  proba-  PRB? 

bility. 


BRG=est , 


RNG=est . 


DEF  A 


ENTER 


ENTER 


ENTER 


ENTER 


DEF  Z 


ENTER 


DEF  X 


ENTER 
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Computed  value  of  the  ellipse  size: 
See  Note  2. 


SIZE=val. 


ENTER 


Step 

Instruction 

Prompt 

Press 

13 

Computed  semi -major  axis  length: 

SMJ-val. 

ENTER 

14 

Computed  major  axis  direction: 

DRC=val. 

ENTER 

15 

Computed  semi-minor  axis  length: 

SMI=val. 

ENTER 

16 

Computed  containment  ellipse  area: 

(See  Note  3.) 

A  =  val. 

17 

For  a  containment  ellipse  of  a  given 
size,  press  DEF,  S. 

DEF  S 

18 

Key  in  the  desired  containment  ellipse 
size  (see  Note  2) . 

SIZE? 

ENTER 

19 

Computed  value  of  the  containment 
probability: 

PRB=val . 

ENTER 

20 

Computed  semi-major  axis  length: 

SMJ=val. 

ENTER 

21 

Computed  major  axis  direction: 

DRC=val . 

ENTER 

22 

Computed  semi-minor  axis  length: 

SMI=val . 

ENTER 

23 

Computed  containment  ellipse  area: 

(See  Note  3.) 

A  =  val. 

24 

To  enter  additional  data  from  new  or 
old  stations,  press  DEF,  C  and  then 
repeat  Steps  2,  3,  4  &  5. 

DEF  C 
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Notes: 


1. 


2. 


3. 

4. 


Reciprocal  bearings  are  used  when  bearings  are  taken  on  known 
positions  from  the  unknown  position  (target) . 

In  the  model  that  is  the  basis  for  the  program,  for  a  given 
probability  of  containment,  the  minimum  area  containment 
region  is  an  ellipse  centered  on  the  position  estimate.  The 
semi-major  axis  =  ka„_  and  semi-minor  axis  =  ko„_  where  k  is 

MJ  Ml 

the  size  of  the  ellipse  and  oMJ  and  oMI  are  the  standard 
deviations  (uncertainty  measure)  of  the  position  estimate 
in  the  major  axis  and  minor  axis  directions. 

The  area  units  are  the  range  units  squared. 

For  a  definition  of  the  program  initiating  keys  and  their 
function,  press  DEF,  H.  The  display  will  show: 

TPE  *  A  EST  =  Z  SIZ  =  S  PRB  *  X.  Next  press  ENTER.  The 
display  will  show:  CON  =  C.  To  repeat  the  display,  press 


ENTER. 

5.  For  a  PC-1211,  use  the  DEF  mode  and  substitute  SHFT  for 
DEF  wherever  DEF  appears  in  the  instructions. 
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Position  Estimation  Program 
FX-702P  User  Instructions 


Step 

Instruction 

Prompt 

Press 

1 

Enter  the  program  (see  Note  1) .  To  run 
the  program,  first  press  FI,  0. 

FI  0 

2 

Next,  key  in  1  and  press  EXE. 

OPTION? 

1  EXE 

3 

Key  in,  in  decimal  degress,  the  obser¬ 
ved  bearing  of  the  target  from  a 
station  or  the  reciprocal  of  the  ob¬ 
served  bearing  of  the  station  from 
the  target  (see  Note  2) . 

OBS  BRG? 

EXE 

4 

Key  in  the  station  bearing  in  decimal 
degrees  from  the  reference  point.  Use 
zero  if  the  station  is  the  reference 
point. 

STA  BRG? 

EXE 

5 

Key  in  the  station  range  (in  any  units) 
from  the  reference  point.  Use  zero 
if  the  station  is  the  reference  point. 

STA  RNG? 

EXE 

6 

Key  in  the  bearing  error  (standard 
deviation)  in  decimal  degrees. 

BRG  ERR? 

EXE 

7 

Repeat  Steps  3 ,  4 ,  5  &  6  for  one  or 
more  additional  stations. 

8 

To  compute  bearing  &  range  estimates, 
first  press  FI,  0. 

FI  0 

9 

Next,  key  in  2  and  press  EXE. 

OPTION? 

2  EXE 

10 

A  computed  target  bearing  estimate  in 
decimal  degrees  is  displayed. 

BRG=est. 

CONT 

11 

A  computed  target  range  estimate  is  dis¬ 
played.  For  an  elliptical  contain¬ 
ment  region,  go  to  Step  12  or  Step  20. 

To  enter  additional  data  from  new  or 
old  stations,'  go  to  Step  28. 

RNG=est. 

12 

For  a  containment  ellipse  of  a  given 
containment  probability,  first  press 

FI/  0. 

FI  0 

13 

I 

Next,  key  in  3  and  press  EXE. 

OPTION? 

3  EXE 

10 


Step 

Instruction 

Prompt 

Press 

14 

Key  in  the  desired  probability. 

PRB? 

EXE 

15 

Computed  value  of  the  ellipse  size: 

See  Note  3. 

SIZE=val. 

CONT 

16 

Computed  semi-major  axis  length: 

SMJ=val. 

CONT 

17 

Computed  major  axis  direction: 

DRC-val. 

CONT 

18 

Computed  semi-minor  axis  length: 

SMI=val . 

CONT 

19 

Computed  containment  ellipse  area: 

See  Note  4 . 

A  =  val. 

20 

For  a  containment  ellipse  of  a  given 
size,  first  press  Fl,  0. 

Fl  0 

21 

Next,  key  in  4  and  press  EXE. 

OPTION? 

4  EXE 

22 

Key  in  the  desired  containment  ellipse 
size  (see  Note  3) . 

SIZE? 

EXE 

23 

Computed  value  of  the  containment 
probability: 

PRB=val. 

CONT 

24 

Computed  semi-major  axis  length: 

SMJ=val . 

CONT 

25 

Computed  major  axis  direction: 

DRC=val . 

CONT 

26 

Computed  semi-minor  axis  length: 

SMI = val. 

CONT 

27 

Computed  containment  ellipse  area: 

See  Note  4 . 

A  =  val. 

28 

To  enter  additional  data  from  new  or 
old  stations,  first  press  Fl,  0. 

Fl  0 

29 

Next,  key  in  5,  press  EXE  and  then 
repeat  Steps  3,  4,  5  &  6. 

OPTION? 

5  EXE 

11 


Notes: 
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POSITION  ESTIMATION  PROGRAM 
TI-59  USER  INSTRUCTIONS 


Step 

Instruction 

Enter 

Press 

Display 

1 

Enter  the  program  (see  Note  1) . 

2 

To  run  the  program,  press  A. 

A 

9 

3 

Enter,  in  decimal  degrees,  the  observed 
bearing  of  the  target  from  a  station  or 
the  reciprocal  of  the  observed  bearing 
of  the  station  from  the  target  (see 

Note  2) . 

e 

B 

0 

4 

Enter  the  station  bearing  in  decimal 

a 

R/S 

a 

degrees  from  the  reference  point.  Use 
zero  if  the  station  is  the  reference 
point . 

5 

Enter  the  station  range  from  the  refer¬ 
ence  point.  Use  zero  if  the  station  is 
the  reference  point.  Use  any  units. 

P 

R/S 

P 

6 

Enter  the  bearing  error  (standard  devi¬ 
ation)  in  decimal  degrees.  After  press¬ 
ing  R/S,  the  display  indicates  the  order 
number  of  the  data  entry. 

e 

R/S 

n 

7 

Repeat  Steps  3,  4,  5  &  6  for  at  least 
one  more  station. 

8 

Display  a  bearing  estimate  $  with  respect 
to  the  reference  station. 

C 

<P 

9 

Display  a  range  estimate  r  with  respect 
to  the  reference  station.  For  an  ellip¬ 
tical  containment  region,  go  to  Step  10 
or  Step  15.  To  enter  additional  data 
from  new  or  old  stations,  repeat  Steps  3, 

4,  5  &  6. 

R/S 

A 

r 

10 

For  a  containment  ellipse  of  a  given  con¬ 
tainment  probability,  enter  the  contain¬ 
ment  probability.  Next,  press  E  and 
display  the  ellipse  size  (see  Note  3) . 

P 

E 

k 

11 

Display  the  semi-major  axis  length. 

R/S 

koMJ 

12 

Display  the  major  axis  direction. 

R/S 

Y 

13 

Display  the  semi -minor  axis  length. 

R/S 

,WMI 

13 


Step 


Instruction 


Enter 


Press 


Display 


. 

14 

Display  the  containment  ellipse  area. 

(See  Note  4.) 

R/S 

area 

15 

For  a  containment  ellipse  of  a  given 
size,  enter  the  containment  ellipse 
size  (See  Note  3) .  Next,  press  D  and 
display  the  containment  probability. 

k 

D 

P 

16 

Display  the  semi-major  axis  length. 

R/S 

k°MJ 

17 

Display  the  major  axis  direction. 

R/S 

y 

18 

Display  the  semi-minor  axis  length. 

R/S 

k°MI 

19 

Display  the  containment  ellipse  area. 

(See  Note  4.) 

R/S 

area 

20 


To  enter  additional  data  from  new  or 
old  stations,  repeat  Steps  3,  4,  5 
&  6  . 


Notes: 


1.  The  program  requires  the  normal  partition.  If  the  calculator 
has  been  in  use,  this  can  be  assured  by  turning  the  calculator 
off  and  then  on  before  loading  the  program. 

2.  Reciprocal  bearings  are  used  when  bearings  are  taken  on  known 
positions  from  the  unknown  position  (target) . 

3.  In  the  model  that  is  the  basis  for  the  program,  for  a  given 
probability  of  containment,  the  minimum  area  containment  region 


is  an  ellipse  centered  on  the  position  estimate.  The  semi¬ 
major  axis  =  ko^  and  semi-minor  axis  =  koMI  w^ere  ^  is  the 
size  of  the  ellipse  and  and  aMI  are  the  standard  deviations 
(uncertainty  measure)  of  the  position  estimate  in  the  major 


axis  and  minor  axis  directions. 


4.  The  area  units  are  the  range  units  squared. 

5.  Negative  bearing  estimates  and  negative  major  axis  directions 
can  result.  To  convert  a  negative  bearing  estimate  to  the 
value  that  would  be  output  by  the  other  programs,  add  360°  to 
the  estimate.  For  example,  -5°  becomes  355°.  To  convert  a 
negative  direction,  add  180°.  For  example,  -5°  becomes  175°. 
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III.  Two  Examples 

In  Scenario  1,  the  scenario  for  the  first  example,  bearings 
are  taken  on  a  target  from  three  separate  stations.  Figure  1  on 
Page  18  shows  Scenario  1  and  Table  1  below  gives  the  station  data. 
The  stations  are  numbered  according  to  the  order  of  station  data 
input  to  the  program. 


TABLE  1 

OBS  BRG 

STA  BRG 

STA  RNG  BRG 

ERR 

Station  1 

038° 

334° 

13500 

4° 

Station  2 

324° 

050° 

11350 

3° 

Station  3 

003° 

000° 

0 

4° 

Note  from  Table 

1  that  the 

reference 

point  is  at  Station  3, 

> 

Program  outputs  for  Scenario  1  are  given  in  List  1  on  Page 
19.  List  1  and  List  2  which  gives  program  outputs  for  Scenario 
2,  are  copies  of  printer  tapes  that  were  generated  with  a  Casio 
FP-10  printer  and  a  Casio  FX-702P  using  the  Casio  FX-702P  program. 
With  allowance  for  differences  in  user  instructions,  display  for¬ 
mats  and  round-off  errors,  the  tapes  indicate  the  output  that 
should  be  obtained  using  any  of  the  other  calculator  programs  ex¬ 
cept  for  the  TI-59  program  where  equivalent  negative  angles  occur. 

In  Scenario  2,  the  scenario  for  the  second  example,  bearings 
are  taken  from  a  target  on  three  separate  stations.  Figure  2  on 
Page  20  shows  Scenario  2  and  Table  2  below  gives  the  station  data. 
As  in  Scenario  1,  the  stations  are  numbered  according  to  the  order 
of  station  data  input  to  the  program.  In  this  scenario  the  reci¬ 
procal  of  the  observed  bearings  are  used  in  order  to  provide  an 
equvalent  scenario  that  is  appropriate  for  the  program. 
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TABLE  2 


OBS  BRG 

RCP  BRG 

STA  BRG 

STA  RNG 

BRG  ERR 

Station  1 

211° 

031° 

000° 

0 

3° 

Station  2 

172° 

352° 

115° 

11100 

3° 

Station  3 

146° 

326° 

082° 

13800 

3° 

Mote  from 

Table  2  that 

the  reference 

point  is 

at  Station 

1. 

Program  outputs 

for  Scenario  2 

are  given 

in  List  2 

on  Page 

21.  The  two  bearing  estimates  in  the  list  are  the  reciprocal  of 
the  target's  estimate  of  the  bearing  of  Station  1. 
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o 


Position  Estimate 


Table  1  gives  the  bearings  of  the  indicated  bearing 
lines  and  List  1  gives  the  bearing  and  the  range  of 
the  position  estimate  from  Station  3. 


LIST  1 


OPTION? 

OPTItN? 

1 

5 

7PE*1  E$T*2  SI 2*3 

DOS  BRG? 

3 

PRfi*4  C0H=5 

STfi  BR6? 

0 

OPTION? 

STftRNG? 

1 

0 

OBSBRG? 

BRG  ERR? 

38 

4 

Sift  BRG? 

OBSBRG? 

334 

$T RUNG? 

OPTION? 

13500 

2 

BOG  ERR? 

BRG*  0.15 

4 

OBSBRG? 

RN6*  19553.70 

324 

STABRG? 

END 

50 

STRONG? 

OPTION’ 

11350 

3 

BOG  ERR? 

SIZE? 

3 

2 

DBS  BRG? 

PRB*  0.86 

OPTION? 

SHJ*  1712.95 

2 

BRG*  359.51 

MR*  12.48 

RH6*  19494.39 

SHI*  1129.98 

END 

fi=  60B0418.40 

OPTION? 

END 

3 

SIZE? 

OPTION? 

2 

4 

PRB=  0.86 

PRB? 

.9 

SHJ*  1137.32 

SIZE*  2.15 

MR*  11.G9 

SHJ*  1837.97 

SHI*  1232.96 

MR*  12.48 

0*  6729444.91 

SHI*  1212.36 

END 

8*  7080340.38 

OPTION? 

END 

4 

PRO? 

.9 

SIZE*  2.15 

SIW*  1064.12 

DIR*  17.69 

SHI*  1322.94 

0*  7747559.77 

END 
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List  2 


OPTION? 

OPTION? 

0 

5 

TPM  £ST=2  SIZs3 

OBS  BRG? 

326 

STfi  BRG? 

82 

STfi  RNfi? 

PRB=4  C0H=5 

OPTION? 

1 

13888 

DBS  BR6? 

BRG  ERR? 

31 

3 

Sift  8R6? 

A 

OBS  BRG? 

STfi  RN6? 

OPTION? 

a 

2 

BRG  ERR? 

J 

OBS  BRG? 

352 

BR6=  31.65 

RN6=  13765.38 

STfi  BRG? 

115 

END 

STfi  RN6? 

OPTION? 

HIM 

3 

BRG  ERR? 

SIZE? 

J 

2 

OBS  BRG? 

PRB=  8.86 

OPTION’ 

2 

BRG=  31.88 

SHJ=  2023.19 

DIR*  172.86 

RN6=  14792.53 

SHI=  967.90 

ENO 

fis  6151983.32 

OPTION? 

3 

SIZE? 

END 

OPTION? 

2 

4 

PRB=  8.86 

PR8? 

.9 

SHJs  3686.56 

SIZEs  2.15 

DIR»  15.39 

SHJs  2178.85 

SHI=  1253.72 

OIRs  172.86 

8=  14205078.91 

$HI=  1838.54 

END 

ft=  7882732.54 

OPTION? 

4 

PRfi? 

.9 

SIZE*  2.15 

SHJs  3869.77 

OIR*  15.39 

$Hls  1345.22 

8=  16354192.26 

END 

END 
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HP-41CV  PROGRAM 


01*LBL  *TPE' 

02  BEG 
03  CLRG 
04  .001 
05  STO  09 
06  26 
07  STO  10 
08*LBL  'CON* 

09  FIX  0 

10  *0BS  BRG?" 

11  PROMPT 

12  STO  11 

13  "STO  BRG?' 

14  PROMPT 

15  STO  12 

16  “STO  RNG?’ 

17  PROMPT 

18  STO  13 

19  'BRG  ERR?" 

20  PROMPT 

21  STO  14 

22  1 

23  ST+  08 

24  2 

25  RCL  08 

26  X>Y? 

27  G TO  02 

28  1 

29  ST+  10 

30  RCL  11 

31  STO  IHD  10 

32  1 

33  ST+  10 

34  RCL  12 

35  STO  IHD  10 

36  1 

37  ST+  10 

38  RCL  13 

39  STO  INB  10 

40  1 

41  ST+  10 

42  RCL  14 

43  STO  IND  10 

44  1 

45  RCL  00 

46  X=Y? 

47  GTO  'COM* 

48  RCL  29 

49  RCL  28 
58  RCL  27 


51  - 

52  SIM 

53  * 

54  STO  15 

55  RCL  33 

56  RCL  32 

57  RCL  31 

58  - 

59  SIN 

60  * 

61  STO  16 

62  RCL  31 

63  RCL  27 

64  - 

65  SIN 

66  STO  17 

67  X=0? 

68  GTO  08 

69  RCL  15 
78  RCL  31 

71  SIN 

72  * 

73  RCL  16 

74  RCL  27 

75  SIN 

76  * 

77  - 

78  RCL  17 

79  / 

88  STO  01 

81  RCL  15 

82  RCL  31 

83  COS 

• 

85  RCL  16 

86  RCL  27 

87  COS 

88  « 

89  - 

98  RCL  17 

91  / 

92  STO  02 

93  26 

94  STO  10 
95»LBL  05 

96  1 

97  ST+  10 

98  RCL  IHD  10 

99  STO  11 

100  1 


101  ST*  10 

102  RCL  INB  10 

103  STO  12 

104  1 

105  ST*  10 

106  RCL  INB  10 

107  STO  13 

108  1 

189  ST*  10 

110  RCL  INB  10 

111  STO  14 

112  XEQ  07 

113  ISG  09 

114  GTO  05 

115  GTO  'CON* 

1 16*LBL  02 

117  XEO  07 

118  GTO  'CON' 

1 19*LBL  07 

120  RCL  01 

121  RCL  13 

122  RCL  12 

123  SIM 

124  * 

125  - 

126  RCL  02 

127  RCL  13 

128  RCL  12 

129  COS 
138  * 

131  - 

132  R-P 

133  STO  16 

134  XOY 

135  STO  15 

136  RCL  11 

137  RCL  15 

138  - 

139  180 
148  XOY 

141  XOY? 

142  GTO  11 

143  360 

144  - 

145*LBL  11 

146  RCL  14 

147  / 

148  STO  17 

149  RCL  16 

150  RCL  14 


23 


f 


151  « 

281  Ra  85 

251  COS 

152  PI 

282  RCL  86 

252  Xt2 

IS  * 

283  * 

253  RCL  85 

154  188 

284  - 

254  * 

IS  / 

285  Ra  28 

255  RCL  23 

IS  STO  18 

286  / 

256  1 

is  «a  15 

287  ♦ 

257  P-R 

IS  COS 

288  STO  26 

258  * 

155  RCL  18 

289  Ra  82 

259  2 

168  / 

218  Ra  83 

268  * 

181  STO  19 

211  RCL  87 

261  RCL  84 

182  RCL  15 

212  * 

262  * 

IS  SIN 

213  RCL  84 

263  STO  26 

164  RCL  18 

214  RCL  86 

264  - 

IS  / 

215  * 

265  RCL  23 

168  STO  28 

216  - 

266  SIH 

167  Ra  19 

217  RCL  28 

267  Xt2 

IS  Xt2 

218  / 

268  RCL  83 

169  ST*  83 

219  * 

269  * 

178  RCL  19 

228  RCL  26 

278  ♦ 

171  RCL  28 

221  XOY 

271  RCL  28 

172  « 

222  R-P 

272  CHS 

173  ST*  84 

223  STO  21 

273  / 

174  Ra  28 

224  XOY 

274  STO  24 

175  Xt2 

225  X>8? 

275  RCL  23 

176  ST*  85 

226  CTO  86 

276  SIH 

177  Ra  17 

227  368 

277  Xt2 

178  Ra  19 

228  ♦ 

278  RCL  85 

179  * 

229*LBL  86 

279  * 

188  ST*  86 

238  STO  22 

288  RCL  26 

181  RCL  17 

231  Ra  84 

281  * 

182  Ra  28 

232  Ra  83 

282  RCL  23 

183  * 

233  RCL  85 

283  COS 

184  ST*  87 

234  X*Y? 

284  XT2 

185  RTH 

235  GTO  83 

285  RCL  83 

t86*LBL  'CST* 

236  RCL  84 

286  ♦ 

187  FIX  2 

237  SIGH 

287  + 

188  RCL  84 

238  45 

288  RCL  28 

189  Xt2 

239  ♦ 

289  CHS 

198  RCL  83 

248  GTO  84 

298  / 

191  Ra  85 

241 *LBL  83 

291  STO  25 

192  * 

242  - 

292  RCL  24 

193  - 

243  / 

293  XOY? 

194  STO  28 

244  2 

294  GTO  81 

195  X=8? 

245  * 

295  STO  25 

196  CTO  88 

246  OTftH 

296  XOY 

197  RCL  81 

247  2 

297  STO  24 

198  RCL  84 

248  / 

298  98 

199  Ra  87 

249*LBL  84 

299  ST*  23 

288  * 

258  STO  23 

388*LBL  81 

24 


381  *BRG=* 

351  * 

382  ARCL  22 

352  RCL  11 

383  AV1EH 

353  RCL  24 

384  STOP 

354  SORT 

385  “RHG=* 

355  * 

386  ARCL  21 

356  XOY 

387  flVIEH 

357  ‘SHJ*’ 

388  STOP 

358  ARCL  X 

388  ’EM’ 

359  flVIEH 

318  GTO  88 

368  STOP 

311*181  -SI2- 

361  RCL  23 

312  'SIZE* 

362  X>8? 

313  PROMPT 

363  CTO  18 

314  STO  11 

364  188 

315  Xt2 

365  * 

316  2 

366«LBL  18 

317  / 

367  ’BRC=‘ 

318  CHS 

368  flRCL  X 

319  EtX 

369  flVIEH 

328  1 

378  STOP 

321  XOY 

371  -SHI=’ 

322  - 

372  RBH 

323  STO  88 

373  flRCL  Y 

324  -PRB=* 

374  flVIEH 

325  ARCL  88 

375  STOP 

326  flVIEH 

376  PI 

327  STOP 

377  * 

328  CTO  89 

378  * 

329*LBL  -PR8- 

379  *«=■ 

338  -PRB?' 

388  flRCL  X 

331  PROMPT 

38i*LBL  88 

332  CHS 

382  flVIEH 

333  1 

383  STOP 

334  + 

384  -EHB- 

335  LH 

385  CTO  88 

336  2 

386*LBL  88 

337  ♦ 

387  flVIEH 

338  CHS 

388  STOP 

339  SORT 

389  -HO  SOL' 

348  STO  11 

398  CTO  88 

341  *SIZE=* 

391  END 

342  flRCL  11 

343  flVIEH 

344  STOP 

345*LBL  89 

346  RCL  23 

347  STO  12 

348  RCL  11 

349  RCL  25 

358  SORT 

25 


SHARP  PC-1500  PROGRAM 


102 "A" : CLEAR  :DIM 
A(7) : DEGREE 
15: "C": INPUT  "08S 
BRG?  ";P: 

INPUT  "STA  BRG 
?  "jQ: INPUT  “S 
TA  RNG?  " >  R: 
INPUT  "BRG  ERR 
?  ";0: PAUSE  " 

It 

20: IF  1=2G0T0  50 
25: l=I+l:A(l-l)=P 
: A< 1+1 )-Q: A( !-*■ 
3)=R: A( l+5)=0: 
IF  1=1G0T0  15 
30:X=A(4)*SIN  (A( 
2)-A(0) ) : Y=A(5 
)*SIN  (A(3)-A( 

1 ) ) : Z=SIN  <AC 1 
)-A(0) ) : IF  2=0 
GOTO  150 

35:U=(X*SIN  AC1)- 
Y*SIN  A (0))/Z: 
U=(X*COS  A(l)- 
Y*COS  A(0))/Z 
40: FOR  n=0TO  1 
45:  P=A(I1) :  Q=A(f1+2 
):R=A(I1+4):0=A 
(M+6) : GOSUB  17 
0: NEXT  M: GOTO 
15 

50: GOSUB  170: GOTO 
15 

60: "Z": PAUSE  "  ": 
F=(B*B-A*C): IF 
F=0GOTO  150 
65:X=U+(B*E-C*0)/ 
F:Y=U-KA*E-B*0 
)/F: GOSUB  200 
70: T=SGN  B*45: IF 
A=CGOTO  80 
75:T=.5*ATN  (2*8/' 
(A-C) ) 

80: G=(C*COS  T*COS 
T-2*B*C0S  T* 
SIN  T+A*S1N  T* 
SIN  T)/-F:G=<TG 
85: H=(C*S1N  T*SIN 
T+2*B*C0S  T* 
SIN  T+A*COS  T* 
COS  T  )/-F :  H=-TH 
: IF  H>=GGOTO  3 
5 

90:Z=H:H=G:G=Z:  T= 
T+90 


95: PRINT  "BR6=";J 
: PRINT  "RNG=" 5 
X: GOTO  145 
100: "H":PR1NT  "  TP 
E=A  EST=Z  SIZ= 
S  PRB=X" 

105: PRINT  "CON=C" : 
GOTO  100 

120: "S": INPUT  "SIZ 
E?  " j  S: 0*1-EXP 
(-S*S/2) 

125: PRINT  USING  "4 
4.44"; "PRB=";0 
: USING  : GOTO  1 
35 

130: "X" : INPUT  "PRB 
?  " ; 0: S=T (-2* 
LN  (1-0)): 
PRINT  "SIZE="; 
S 

135: X=S*G: Y-S*H: N= 
T : PRINT  "SnJ=" 

; Y: IF  N<0LET  N 
=N+188 

140:  PRINT  "DIR=‘';N 

: print  "sm="; 

x:  PRINT  ,'A=";  ri 
*X*Y 

145: PRINT  "END": 
GOTO  145 

150: PRINT  "NO  SOL" 

: GOTO  150 

170: X=U-R*SIN  Q:Y= 
U-R*COS  Q: 
GOSUB  200 

175:  U1=(P-J) :  L=K*0* 
JT/180:G=COS  J/ 
L:H=S1N  J/L: IF 
U>=1 80LET  U=W- 
360: GOTO  185 
180: IF  U<=- 1 80LET 
U=W+360 

185: 11=14/0  :A=G*G+ A: 
B=G*H+B:C=H*H+ 
C:0=U*G+D:E=W* 
H+E: RETURN 
200: K=T (X*X+Y*Y) : 
IF  K=0LET  J=0: 
RETURN 

205:  J=ACS  (Y/-K):  IF 
ASN  (X/KX0LET 
J=360-J 
210: RETURN 
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SHARP  PC-1211  PROGRAM 


10* "A"* CLEAR  * 
DEGREE 

15* “C°* INPUT  "0 
BS  BRG?  "IP* 
INPUT  "STA  B 
RG?  "IQ* 
INPUT  "STA  R 
NG?  “IR* 
INPUT  “BRG  E 
R R?  "10* IF  1 
-2G0T0  130 

75* I»I+1*A< 1+26 
>*P*  A< 1+28)* 
Q*A< 1+30 >*R* 
A<  I+32>*0*  IF 
I*1G0T0  15 

85*X*A<31)*SIN 
<A<29)-AC2?> 
)*Y*A<32>* 
SIN  <A<30)-A 
<28>>*Z*SIN 
<A<28)-A<2?> 
> 

90* IF  Z*OGOTO  3 
25 

95*U»<X*SIN  A<2 
8)-Y*SIN  A<2 
7)VZ  sV-CX* 
COS  A<28)-Y* 
COS  A<27>VZ 
105* FOR  1*170  2* 
P=A<I+26>*Q® 
A(I+28)*R*A< 
I+3G)*0*A<I+ 
32) *  GOSUB  40 
0 

110*»€XT  I* GOTO 
15 

130*  GOSUB  400 
135*  GOTO  15 
140* 

♦O*  IF  F=0 
GOTO  325 
1 50  *  X“U+  <  B*E  -C*D 
VF*Y*V+<A*E 
-B*D>^F: 
GOSUB  500 
160*  T*SGN  B*45* 
IF  ACTHEN  1 
70 

165*  T*.5*ATN  <2* 
&/<A-C>> 
170*G»<C*C0S  T* 
COS  T-2*B* 
COS  T+SIM  T+ 
A*SIH  T+SIN 
T)/--F*G*/G 
175*rt»<C*SIN  T* 
SIN  T+2*B* 
COS  T+SIN  T+ 
A+COS  T*C0S 
T)/'-F*H*TH* 
IF  H>*GG0T0 
135 


180*  Z=H:  H=G*  G=Z* 
T=T+5C 

135*  PRINT  "BRG*" 
}J* PRINT  "RN 
G*" IK* GOTO  3 
20 

200* "H" *PRINT  "T 
FE*A  EST*2  S 
IZ*S  PRB*X" 

205*  FRINT  MCON*C 
"sGOTO  200 

225* "3" *  INPUT  "S 
IZE?  "IS*0*1 
-EXP  <-$*S/2 
> 

230* PRINT  USING 
"##,  #Ti  "PRB 
*"50* USING  * 
GOTO  300 

235* “X"! INPUT  "P 
RB?  "50*S=r< 
-2»LN  0-0)) 
SPRINT  "SIZE 
="5S 

300 *  X=S=*  G  *  Y=  S*H  * 
N=T* PRINT  "S 
MJ=" I Y* IF  N 
OOLET  N*N+1 
80 

315* PRINT  "DRC»" 
5N*PRINT  "SN 
I«"5X*Z*H*X* 
Y* PRINT  "A=" 
*  ^ 

320* PRINT  “END"* 
GOTO  320 

325* PRINT  "NO  SO 
L“*GOTO  325 

400*  X*U-R*S1N  Q* 
y»V-R*COS  q* 
GOSUB  500 *W= 
c'P-J  >  i  L*K*0* 
&/130*G-C0S 
X-'L*  H*  SIN  Jy 
L 

415* IF  U>=130LET 
W*W-360* GOTO 
440 

425* IF  WO-130 
LET  W*W+36Q 

440 *  W*W-'  0  J  A*  G*G+ 
A*B*G*H+B*C* 
H*H+C*D*W+G+ 
D*E*W*H+E* 
RETURN 

500*  K*T<X*X+Y*Y> 
*  IF  K*OLET  J 
■0* RETURN 

510*  M*ASN  <X.-'tO  * 
J*ACS  (Y  -'K.)s 
IF  MCOLET  J* 
360-  J 

515* RETURN 


CASIO  FX-702P  PROGRAM 


5  IMP  "OPTION", J; 

SET  F2:N00E  4 
18  IF  J=0  THEN  175 
15  IF  J=2  THEN  188 
28  IF  J=3  THEN  145 
25  IF  J=4  THEN  158 
38  IF  J=5  THEN  45 
48  VAC  iSET  F2 
45  INP  "OSS  BRfi",P 
s  INP  "STA  BRG", 
Q: IMP  "STA  RNG" 
,R 

58  INP  "BAG  ERR*,0 
s IF  1=2  THEN  95 
55  I«WUKH)«P» 
A(I+l)=Q:A<I+3) 
=R:A(I+5>=0 
68  IF  1=1  THEN  45 
65  X=A(4)*SIN  (A(2 
)-A(8)):Y=A<5)* 
SIN  <A(3)-A(D) 
78  Z=$IN  (A(1)-A(0 
))! IF  Z=8  THEN 
178 

75  U=(X*SIN  A(l)-Y 
♦SIN  A(0))/Z:V= 
<X*COS  A( 1 )-Y*C 
OS  A<8)  )'Z 
80  FOR  H=0  TO  l:P= 
A(H):8=A(H+2) 

85  R=A(«+4)JO=A(N4 
6):GSB  200! NEXT 
Ns  GOTO  45 
90  R=A<«+4)!0=A(H+ 
6):6SB  200! NEXT 
HiGOTO  45 
95  GSfi  200! GOTO  45 
110  F=(0*0-A*C):IF 
F*8  THEN  178 
185  X*U+<B*E-CtO)/F 
:Y=V+<A*£-B*D)/ 
FsRPC  Y.X 

110  T=SGN  B*45! IF  A 
=C  THEN  120 
115  T=.5*RTN  (2*8/< 
(K)) 

120  G*(C*COS  T*COS 
T-2*8*C0S  T*SIN 
T+A*SIN  T*SIH 
T>/-F 


125  H»<C*SIN  T*SIN 
T+2*fi*COS  TtSIN 
T+A*CQS  T*COS 
T)/-F 

138  6*S6R  G!H=SOR  N 
!IF  B>H;Z=H!H=G 
!G*Z!T=T+90 
135  IF  Y<0JY=Y+368 
148  PRT  ”BR6=*;YiPR 
T  *RNG="!X!SOTO 
165 

145  INP  "SI2E",$:0= 
1-EXP  (-S*S/2): 
PRT  "PRB=*;0:G0 
TO  155 

158  INP  "PRB”,0:S=$ 
OR  (-2*LH  U-Q> 
)iPRT  "SI ZE=";S 
155  X=S*6!Y=S*H:PRT 
”SHJ=";Y:N=T:I 
F  N<0;N=N+I80 
168  PRT  "DIR=*:H:PR 
T  "SHI=»;X:PRT 
"A=';iWY 
165  PRT  "END": GOTO 
165 

170  PRT  "NO  SOL":  GO 
TO  178 

175  PRT  *TPE=1  EST= 
2  SIZ=3” 

188  PRT  "PRB=4  CON= 
5”: GOTO  175 
208  X=U-R*SIN  Q:Y=V 
-R*COS  0:RPC  Y, 
X 

285  «=<P-Y)!L=X*0** 
'180:8=005  Y/L: 
H=SIN  Y/L 

218  IF  Hil80SH=H-36 
8:60T0  228 
215  IF  Hi-188; H=H+3 
68 

228  H=H/Q:A=6*6+A;B 
*G*H+B:C=«*H+C: 
D=H*6+D:E=H*H+E 
:RET 


28 


TI-59  PROGRAM 


000 

47 

CMS 

050 

10 

E* 

100 

42 

STD 

001 

60 

DEG 

051 

18 

C* 

101 

38 

002 

02 

2 

052 

87 

IFF 

102 

43 

RCL 

003 

42 

STD 

053 

03 

03 

103 

14 

14 

004 

09 

09 

054 

01 

01 

104 

39 

CDS 

005 

69 

DP 

055 

28 

28 

105 

65 

X 

006 

28 

28 

056 

19 

D* 

106 

43 

RCL 

007 

76 

LBL 

057 

75 

- 

107 

18 

18 

008 

16 

fl* 

058 

19 

D* 

108 

75 

- 

009 

09 

Q 

059 

95 

= 

109 

43 

RCL 

010 

42 

STD 

060 

94 

+  /- 

110 

10 

10 

Oil 

00 

00 

061 

38 

SIN 

111 

39 

CDS 

012 

42 

STO 

062 

65 

x 

112 

65 

v 

013 

01 

01 

063 

19 

It* 

113 

43 

RCL 

014 

92 

RTN 

064 

95 

= 

114 

19 

19 

015 

76 

LBL 

065 

48 

EMC 

115 

95 

= 

016 

18 

C* 

066 

19 

19 

116 

55 

-T 

017 

69 

□P 

067 

22 

INV 

117 

43 

RCL 

018 

20 

20 

068 

97 

DSZ 

118 

20 

20 

019 

-?  o 
r  £ 

ST* 

069 

09 

09 

119 

y 

SIN 

020 

00 

00 

070 

00 

00 

120 

= 

021 

92 

RTN 

071 

76 

76 

121 

42 

STD 

022 

76 

LBL 

072 

01 

1 

122 

39 

39 

023 

19 

D' 

073 

69 

DP 

123 

86 

sir 

024 

69 

□P 

074 

21 

21 

124 

03 

03 

025 

21 

21 

075 

92 

RTN 

125 

02 

026 

73 

RC* 

076 

42 

STD 

126 

42 

STD 

027 

01 

01 

077 

18 

18 

127 

09 

09 

028 

•92 

RTN 

078 

65 

y 

128 

16 

FT 

029 

76 

LBL 

079 

43 

RCL 

129 

19 

Ii 1 

030 

10 

E' 

080 

14 

14 

130 

10 

E  ' 

031 

65 

X 

081 

42 

STD 

131 

42 

STD 

032 

89 

082 

20 

20 

132 

18 

18 

033 

55 

-r 

083 

38 

SIN 

133 

19 

Ii' 

034 

01 

1 

084 

75 

- 

134 

■— 1  aL- 

XI T 

035 

08 

8 

085 

43 

RCL 

135 

19 

D' 

036 

00 

0 

086 

10 

10 

136 

32 

XI T 

037 

95 

= 

087 

•y 

INV 

137 

37 

P/R 

038 

92 

RTN 

088 

44 

SUM 

138 

75 

- 

039 

76 

LBL 

089 

20 

20 

139 

43 

RCL 

040 

11 

ft 

090 

38 

SIN 

140 

38 

38 

041 

81 

RST 

091 

65 

x 

141 

95 

— 

042 

76 

LBL 

092 

43 

RCL 

142 

94 

+  /  - 

043 

12 

B 

093 

19 

19 

143 

32 

XI T 

044 

18 

C 

094 

95 

= 

144 

75 

- 

045 

91 

R/S 

095 

55 

■f 

145 

43 

RCL 

046 

18 

C 

096 

43 

RCL 

146 

39 

39 

047 

91 

R/S 

097 

20 

20 

147 

95 

= 

048 

18 

C' 

098 

38 

SIN 

148 

94 

+  /  - 

049 

91 

R/S 

099 

95 

= 

149 

32 

V  *  T 

A  <  1 

29 


i 


30 


450 

32 

KIT 

451 

65 

X 

452 

89 

it' 

453 

95 

= 

454 

92 

RTN 

455 

61 

GTD 

45b 

04 

04 

457 

54 

54 

458 

76 

LBL 

458 

14 

D 

480 

33 

V  -• 

461 

55 

462 

32 

X 1 T 

463 

02 

Cm 

464 

95 

465 

94 

+  /  " 

466 

■*> 

IHV 

467 

23 

LNX 

4  6  8 

75 

- 

469 

01 

1 

470 

95 

= 

471 

94 

+ /  - 

472 

91 

R.-S 

473 

32 

XI T 

474 

34 

rx 

475 

61 

GTD 

476 

04 

04 

477 

10 

10 

478 

00 

0 

479 

00 

0 

} 


32 


Appendix  1.  The  Estimation  Procedure 

In  this  development  the  assumptions  stated  in  Section  I  are 
required  conditions.  A  rectangular  coordinate  system  is  used 
with  the  positive  y-axis  directed  north,  the  positive  x-axis 
directed  east  and  the  origin  at  the  reference  point  and  all  angles 
are  in  radians.  Figure  3  shows  a  station  located  with  respect  to 
the  coordinate  system.  The  bearing  line  of  length  r  goes  to  the 
object's  unknown  position,  the  bearing  line  of  length  r  goes  to 
an  initial  estimate  of  the  object's  position  and  the  third  bear¬ 
ing  line  corresponds  to  an  observed  bearing. 

To  determine  the  coordinates  for  a  final  estimate,  consider 
the  arc  segments  u  =  r  ( 0 — <J) )  where  0-<f>  is  the  bearing  error  and 
v  =  r(<(>-6)  and  w  =  r(0-0)  that  are  defined  by  the  three  bearing 
lines  and  the  circle  of  radius  r  that  is  centered  on  the  station 
and  goes  through  the  initial  estimate.  The  geometry  is  shown  in 
Figure  3.  Note  that  u  can  be  defined  by  u  =  w  -  v.  In  this 
expression,  w  =  r(6-0)  is  known,  and  v  can  be  determined  in  terms 
of  x  and  y  the  unknown  coordinates  of  the  target.  With  the  refer¬ 
ence  point  at  the  initial  estimte,  to  first  order, 
v  =  x  cos  0  -  y  sin  0  and  u  =  w  -  x  cos  6  +  y  sin  0.  Since  this 
approximation  applies  to  all  stations,  its  use  suggests  that, 
for  each  station,  r  ~  r  which  is  equivalent  to  having  the  initial 
estimate  relatively  close  to  the  target's  position. 

Since,  for  each  station  i,  an  observed  bearing  0^  is  the 
value  of  a  normal  random  variable  0^  with  mean  ,  the  coordinate 
u^  =  r^(9^-$^)  is  the  value  of  a  normal  random  variable  with 

mean  zero.  In  addition,  since  the  0^_  are  independent,  the  Ui  are 
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North 

I  y 


Initial 

Estimate 


Target 


Figure  3 


.  The  coordinate  geometry.  The  coordinates  of  the  initial 
estimate  are  (x*,y*).  In  the  development,  the  reference 
point  is  at  the  initial  estimate. 


also  independent.  And,  consequently,  the  joint  distribution  of 
the  IK  is  determined. 

To  estimate  x  and  y,  maximum  likelihood  estimates  are  used 
here.  The  likelihood  for  a  sample  8^,  9n  from  n  stations 

is 


L(8^, • • • »  ®n)  — 


n 

n 

1 


1 

/2tT 


exp  - 


1 

2 


n  2  2 
2  (6  .  —  <f> . )  /ef 

3.  1  1  1 


and  the  likelihood  for  a  cooresponding  sample  u^,  un  is 

n  1  1  n  2  2 

L  (u.  ,u-,...,u  )  =  II  -  exp  -  -k  E  u./o. 

12  n  1  /2?  o±  2  1  1  1 

with  ck  =  r^e^  where  is  the  standard  deviation  of  0^. 

By  definition,  the  maximum  likelihood  estimates  for  x  and  y 
are  the  estimates  which  make  LCu-^,  u2,...,un)  a  maximum.  In 
this  case,  making  L {u1,u2, . . . ,un)  a  maximum  is  equivalent  to 

n  2  2 

making  E(u.  /a.  )  a  minimum.  So,  to  find  the  maximum  likelihood 
1  /  \  ^ 
estimates  x  and  y,  solve  the  following  two  equations  for  x  and  y: 


9 ( In  L)  1 
9x  1 

x 


0  and 


A 


X 


9 (in  L)  , 
9y  1 

v 


0 

y 


With  =  wi  -  x  cos  8.^  +  y  sin  0^  the  equations  can  be  written 
as  follows: 


and 


n  A  2 
E  (w^  -  x  cos  6i  +  y  sin  6.^)  (cos  8.^/ck  =  0 

n  2 

E  (w.  -  x  cos  0.  +  y  sin  0.) (sin  0.)/a.  =  0  . 

_  1  X  X  XX 
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In  terms  of  the  following  quantities: 

A  =  E(cos2  0^/a2  ,  B  =  E(sin  6i  cos  0^/a2  , 

C  =  E(sin2  0i)/a2  ,  D  =  E  (w^  cos  B^/a2  , 

E  =  E (wi  sin  6i) /a2  , 

the  equations  become : 

Ax  -  By  =  D 

Bx  -  Cy  =  E 


The  solutions  are: 

(1)  x  =  (BE  -  CD) / (B2  -  AC) 
and 

(2)  y  =  (AE  -  BD)/(B2  -  AC) 

A  confidence  region  can  be  constructed  about  an  estimated 
position.  In  order  to  indicate  how  this  can  be  done,  a  proba¬ 
bility  region  about  the  true  position  will  be  considered  first. 

Both  x  and  y  are  values  of  random  variables.  If  a  new  set 
of  bearings  0-^,  02i  . .. ,  9n  is  observed  (for  the  same  initial 
estimate  and  a  fixed  target) ,  in  general,  a  new  pair  of  values 
x  and  y  will  be  obtained. 
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If  X  and  Y  represent  these  random  variables, 

n  2 

£  (Wi/ai)  (B  sin  Bj.  -  C  cos  8i) 

n  2 

£  (VT/ck)  (a  sin  8i  -  B  cos  8i) 


where  -  6i)  . 

/S  A 

Since  X  and  Y  are  a  linear  combination  of  the  n  normal  random 
variables  W^,  W2,  ...»  Wn,  or  equivalently  of  the  n  normal  random 
variables  0^,  02,  ...»  0n,  they  have  a  joint  normal  distribution. 
Since  E(W^)  =  -  8i)  ,  if  8^  =  for  i  =  1,  2,  ...»  n,  that 

is,  if  the  initial  estimate  of  the  target's  position  is  at  the 
target's  position,  E(W1)  =  0  for  i  =  1,  2,  ...»  n.  In  this  case 
E (X)  =0  and  E(Y)  =0  and  the  joint  normal  distribution  is  centered 
on  the  object's  position.  To  the  degree  of  the  approximations 
that  have  been  made,  this  is  also  true  if  the  initial  estimate  is 
not  at  the  target's  position. 

A  region  of  minimum  area  for  a  given  probability  of  contain¬ 
ment  of  an  estimated  position  can  be  determined.  The  region  is 
bounded  by  an  ellipse  which  is  centered  on  the  object's  position 
and  whose  axes  lie  along  the  axes  of  an  x'y' -coordinate  system 
that  is  obtained  by  rotating  the  xy-coordinate  system  that  is 
centered  on  the  object's  position  through  an  angle  y.  In  this 
system,  a~,  is  0,  that  is  x'  and  y'  are  independent  normal  ran- 
dom  variables.  The  two  coordinate  systems  are  illustrated  in 


x  =  -4 


( B  -AC ) 


Y  = 


(B2-AC) 
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Figure  4.  The  coordinates  of  a  point  in  the  two  systems  are 
related  by 


x' 

=  X  cos 

Y 

-  y 

sin 

Y 

y’ 

=  x  sin 

Y 

+  y 

cos 

Y 

These  relations  imply: 


2  2  2  2  2 

(3)  cu  =  a „  dos  y  -  2o^  cos  y  sin  y  +  aA  sin  y  , 

x '  x  xy  y 


2  2  2  2  2 

(4)  oa  =  0„  sin  y  +  2a^„  cos  y  sin  y  +  cos  y 

y '  x  xy  y 


and 


2  2  2  2  2 

(5)  A  =  (o„  -  oj  sin  y  cos  y  +  0^  (cos  y  -  sin  y) 

x'y'  x  y  xy 


where  y,  the  angle  of  rotation  of  the  coordinate  axes,  is  posi¬ 
tive  in  the  clockwise  direction.  And  o^,^,  =  0  implies 


20AA 

tan  2y  =  -XX 

**  - 
y  x 

With  the  initial  estimate  of  the  target's  position  at  the 

2 

target's  position  E(W^)  =  0  and  therefore,  Var(W^)  =  0^  for 
i  =  1,  2,  ...»  n.  In  this  case 

2  1  n  2  2 

oA  *  — 5 - 5-  E  (l/0.j)  (B  sin  6.  -  C  cos  $.)  , 

x  (B  -AC)  ^  1  1  1  1 
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— 5 - 7 

<B  -AC) 


n 

I 

1 


(1/of)  (A 


sin 


6i  ~ 


B  cos 


3i> 


and 


l  “  2 

oAA  *  — 5 -  £  (l/o.)  (B  sin  3.  -  C  cos  B-)  (A  sin  3-  -  B  cos  3,-)  . 

xy  (B  -AC)  1  1  1  11  1 


Using  the  definition  for  A,  B  and  C,  the  above  become 


(6) 

(7) 


and 


_2  _  C 
0~  “  T~  ' 

X  (AC-B  j 

-2  _  A 
°a  ~  5~  ' 

y  (AC— B  ; 


(8) 


B 

- 5-  • 

(AC-B  ) 


So  tan  2  =  2B/(A-C)  for  3^  =  <f>^,  i  =  1,  2,  ...»  n. 

With  the  target's  position  known  and,  consequently,  <|k  known 

2  2 

for  i  =  1,  2,  ...»  n,  the  above  expressions  for  o a  ,  a a  ,Oaa  and  y 

x  y  xy 

can  be  used,  since  the  initial  estimate  of  the  target's  position 

can  be  taken  as  the  target's  position. 

With  values  for  ,  o^a  and  Y,  values  for  Oa,  and  Oa, 

X  y  xy  x  y 

can  be  found  by  using  equations  (3)  and  (4) .  The  probability 

that  an  estimated  position  will  be  within  an  ellipse  of  semiaxes 

kcr£i  and  ko^,  which  is  centered  on  the  target's  position  is 
x  y 
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2 

1  -exp(-k  /2) .  This  result  can  be  found  by  integrating  the  bi¬ 
variate  normal  density  over  the  ellipse.  And  the  area  of  the 
2 

ellipse  is  k o%oa,. 

x  y 

Given  estimates  x  and  y  found  by  using  Equations  (1)  and  (2) , 

the  ellipse  with  semi-axes  ka^,  and  ka~,  in  a  coordinate  system 

x  y 

that  is  centered  on  the  point  (x,y)  and  has  been  rotated  through 

2 

an  angle  y  is  a  1  -  exp(-k  /2)  confidence  region.  This  follows, 
since,  to  the  degree  of  the  approximations  involved,  the  bivar¬ 
iate  normal  distribution  of  X  and  Y  is  centered  on  the  target's 
position.  The  confidence  ellipse  is  defined  if  and  c? ~ ~ 

can  be  found,  that  is  if  the  elements  of  the  covariance  matrix 
can  be  found.  To  the  degree  of  the  approximations  involved, 
this  can  be  done  as  follows:  First,  assume  the  initial  estimate 
of  the  target's  position  is  at  the  target's  position.  Then, 
values  for  cr~  ^  a*  f  and  y  can  be  determined  by  using  Equations 

(6),  (7)  and  (8).  These  values  can  then  be  used  to  determine 
2  2 

°x'  '  °y'  and  °x'y'  using  Equations  (3),  (4)  and  (5).  Now, 
with  a  value  for  k,  a  confidence  region  can  be  constructed.  To 
the  degree  of  the  approximations  involved,  the  shape  of  the  con¬ 
fidence  region  is  independent  of  both  the  target's  position  and 
of  the  initial  estimate  of  the  target's  position. 

For  the  case  where  bearings  are  taken  from  the  target  on  two 
or  more  stations,  0^  is  the  reciprocal  of  the  bearing  taken  from 
the  target. 

A  discussion  of  the  theory  of  bearings  only  position 
estimation  procedures  for  situations  similar  to  the  one  considered 
here  is  given  in  Reference  1.  Reference  2  gives  an  essentially 
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equivalent  bearings  only  procedure.  It  also  gives  a  range  only 
procedure,  a  range  and  bearing  procedure  and  HP-9830A  programs 
with  which  to  implement  the  procedures.  Using  the  fix  determined 
by  two  lines  of  bearing  as  the  initial  estimate  was  suggested  by 
this  reference. 

The  following  equations  are  evaluated  in  the  program  to  deter¬ 
mine  the  coordinates  of  the  initial  estimate: 

x*  sin  (82“^)  =  [p^  sin  {ct^-0^)]  sin  02 

-  [p2  sin  (a2~e2)]  sin  0.^ 

and 

y*  sin  =  ^pl  s*n  ^  003  ®2 

-  [p2  sin  (a2-02)]  cos  0^  . 

Reference  3  describes  a  Tl-59  program  that  is  based  on  an  equi¬ 
valent  procedure.  The  program  allows  a  user  to  either  input 
the  coordinates  of  the  initial  estimate  or  determine  them  in 
the  manner  described  here. 
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Appendix  2.  HP-41CV  Program  Labels 

The  global  labels  in  the  HP-41CV  program  that  are  assigned 
to  the  keys  I+,  1/x,  /x,  LOG  and  LN  give  these  keys  a  mnemonic 
character.  For  example,  with  the  calculator  in  USER  mode,  if 
LOG  is  pressed  and  held,  SIZ  will  be  displayed  and  then  after 
a  delay,  NULL. 

If  a  global  label  is  replaced  by  a  local  label,  this  mnemonic 
character  will  be  lost.  However,  automatic  key  assignment  will 
be  gained  if  the  key  corresponding  to  the  local  label  has  not 
been  previously  assigned.  If  this  is  the  case  and  if  there  is 
a  second  program  in  program  memory  that  uses  the  same  local  label, 
then  it  will  be  automatically  assigned  when  the  calculator  is 
positioned  in  program  memory  at  that  program. 

Global  labels  in  the  HP-41CV  program  can  be  replaced  with 
local  labels.  If  this  is  done  as  described  below,  the  user  in¬ 
structions  will  still  be  applicable.  First,  make  the  following 
replacements : 


Line 

329: 

LBL 

"PRB  " 

LBL 

E. 

Line 

311:  LBL  "SIZ"  -*• 

LBL  D. 

Line 

186: 

LBL 

«est" 

-► 

LBL 

B. 

Line 

008:  LBL  "CON"  -*• 

LBL  C. 

Line 

047: 

GTO 

"CON" 

GTO 

C. 

Next, 

after 

line 

001: 

LBL 

"TPE" 

,  insert  LBL  A  so 

that  it 

becomes  line  002.  Finally,  if  there  is  a  key  assignment  for  a 
key  in  the  following  list:  E+,  1/x,  /x,  LOG  and  LN,  remove  it. 
Now,  with  the  calculator  in  USER  mode,  if  LOG  is  pressed  and  held, 
XEQ  D  will  be  displayed  and  then,  after  a  delay,  NULL.  The  key's 
mnemonic  character  has  been  lost,  but  if  there  is  a  second  pro¬ 
gram  in  program  memory  with  the  local  label  D,  LOG  will  be  auto¬ 
matically  assigned  to  D  when  the  calculator  is  positioned  in  pro¬ 
gram  memory  at  that  program. 
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